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Abstract 


A  system  has  been  developed  which  allows  the  photo  interpreter  to  determine 
the  tilt  of  an  aerial  photograph.  Input  consists  of  image  measurements  made 
using  a  calibrated  template  and  tube  magnifier,  and  three  geodetic  control 
points  scaied  from  large  scale  map  source.  Church's  resection  is  used  to 
compute  the  camera  parameters  on  the  Hewlett-Packard  HP  <HC  program¬ 
mable  calculator.  Tilt  determination  allows  computation  of  lengths,  widths, 
and  heights  from  a  single  image.  These  calculations  follow  naturally  from  the 
resection  and  are  included  in  the  system. 

The  system  incorporates  several  desirable  features.  It  is  easy  to  learn  and  to 
use,  and  the  resection  converges  within  five  minutes.  The  basic  equipment  is 
inexpensive  and  the  calculator  may  be  used  for  other  applications. 


i.  INTRODUCTION: 


A  photo  interpreter  or  photogram metrist  occasionally  is  presented  non¬ 
vertical  or  oblique  photography ;  as  his  primary  source  for  photogrammetric 
measurement.  Poor  flying  conditions  or  inexpert  photographic  technique  fre¬ 
quently  result  in  undesired  nonverticality.  Oblique  photography  may  be  chosen 
because  a  larger  area  of  coverage  can  be  achieved,  and  a  look  at  the  side  of  a 
feature  to  be  analyzed  permits  a  more  reliable  interpretation. 

Whatever  the  source  of  tilt,  the  relatively  simple  measurements  of  vertical 
photography  become  extremely  complicated  and,  without  the  assistance  of  a 
computer,  prohibitively  time-consuming. 

A  variety  of  disciplines  rely  upon  on-the-site  field  measurements.  For 
many,  the  level  of  sophistication  does  not  warrant  full-time  use  of  a  non¬ 
portable  computing  system.  Recent  developments  in  handheld  programmable 
calculators  make  it  possible  to  quickly  resolve  the  tedious  computations  needed 
to  determine  the  tilt  of  an  aerial  photograph.  These  instruments  are  precise 
and  easy  to  use.  Simple  operating  instructions  and  programmed  prompts  lead 
the  operator  through  the  procedure  step-by-step  with  virtually  no  opportunity 
for  blunders. 

A  computer  program  (Appendix  B)  was  developed  at  the  Defense  Mapping 
Agency  Aerospace  Center  (DMA AC)  for  the  Hewlett  Packard  41C  calculator 
which  determines  the  tilt  of  an  aerial  photograph  using  Church's  Resection 
(Church,  1945).  Horizontal  and  vertical  ground  measurements  may  then  be 
computed  for  features  on  that  photograph. 

2.  SOFTWARE: 

Church's  Resection  was  chosen  because  the  math  model,  computation 
algorithm,  and  pre-tested  data  were  readily  available.  Professor  Earl  Church 
developed  the  procedure  and  published  it  in  1945  in  a  Syracuse  University 
technical  paper.  Virtually  every  theoretical  photogrammetry  student  has 
cranked  out  the  solution  on  a  desk  calculator  or  has  programmed  it  in 
FORTRAN.  The  desk  calculator  procedure  reputedly  takes  about  six  hours, 


assuming  no  blunders.  Present  day  students,  after  struggling  through  the  first 
iteration,  opt  for  the  FORTRAN  version  (possibly  an  exercise  designed  to  make 
students  more  receptive  to  the  collinearity  and  coplanarity  equations). 

Professor  Church's  solution  is  iterative,  and  uses  matrices  and  determi¬ 
nates  to  derive  the  exposure  station  coordinates  and  the  swing  and  tilt  of  the 
aerial  photograph  (see  American  Society  of  Photogrammetry  Manual  of  Photo^ 
grammetry,  1980).  These  values  are  defined  in  a  rather  unique  coordinate 
system.  He  presents  the  solution  in  a  "cookbook"  form,  and  provides  checks 
throughout  to  identify  blunders. 

Heights  (Appendix  A.II)  and  lengths  (Appendix  A.UI)  are  computed  by  three 
digital  rectification  algorithms.  One  subroutine  computes  the  direction  cosines 
(Appendix  A.I)  of  vectors  from  the  principal  point  to  the  imaged  nadir  and  to 
the  base  and  top  or  the  end  points  of  the  feature  being  measured.  Image 
coordinates  are  then  projected  onto  a  horizontal  plane  passing  through  the  nadir 
image.  The  vector  passing  from  the  principal  point  through  the  nadir  is 
substituted  for  focal  length  in  the  conventional  single  image  height  and  length 
equations  (figures  1,  2,  and  3). 

The  "software  design  features  repeated  use  of  the  direction  cosine 
algorithm  to  reduce  programming  steps  and  conserve  storage. 

3.  HARDWARE: 

The  HP  41C  programmable  calculator  (Hewlett-Packard,  1979)  was  chosen 
because  it  had  memory  expandable  to  the  requirements  of  the  solution  and 
because  of  the  authors'  familiarity  with  the  programming  language.  (Choice  of 
this  calculator  does  not  constitute  endorsement  by  DMAAC  of  Hewlett  Packard 
or  the  HP  41C  calculator.)  The  HP  41C  features  continuous  memory,  LED 
display,  and  magnetic  card  reader.  With  memory  modules,  the  calculator  is 
expandable  to  more  than  300  registers  (250  with  card  reader  attached).  It  is 
battery  operated  with  very  low  power  consumption,  and  accessories  include 
rechargeable  batteries.  The  system  is  suitable  for  office  or  field  applications. 


4.  PROCEDURES 


Computing  heights  or  lengths  with  this  system  consists  of  measuring  map 
coordinates,  measuring  image  .coordinates,  and  operating  the  calculator 
program.  The  system  concept  was  developed  to  increase  accuracies  with  a 
minimum  of  operator  training  and  a  minimum  of  opportunity  for  operator  error. 

Map  derived  control  is  scaled  from  large  scale  map  source  such  as  United 
States  Geological  Survey  (USGS)  1:24000  scale  topographic  maps.  The  user 
reads  the  Northing  and  Easting  from  the  Universal  Transverse  Mercator  (UTM) 
grid  with,  for  example,  the  Metric  Coordinate  Reader  (CR-2)  (USGS,  1973), 
manufactured  by  USGS.  The  Geological  Survey  reports  that  15  meter  accu¬ 
racies  are  expected  with  their  1:24000  maps. 

A  1-centimeter  grid  aligned  to  the  fiducials  of  the  photo  is  fixed  in  place. 
To  measure  photo  coordinates,  the  user  counts  the  grid  lines  to  determine 
centimeter  values.  Usmg  a  Bausch  &  Lomb  calibrated  tube  magnifier,  he  reads 
x  and  y  coordinates  to  a  least  count  of  0.1  millimeter.  (Choice  of  the  Bausch  & 
Lomb  calibrated  tube  magnifier  does  not  constitute  endorsement  by  DMAAC  of 
Bausch  &  Lomb  or  its  products.)  Papa-  stretch  of  the  photo,  image  quality,  and 
other  factors  will  contribute  to  error;  however,  0.2  millimeter  accuracies  are 
expected  as  a  worst  case. 

To  operate  the  calculator  program,  the  user  first  runs  the  routines  for 
Church's  Resection  (Appendix  B,  Line  01).  Algorithms  for  length  or  height 
determination  are  then  applied. 

Church's  Resection  requires  as  input:  1)  UTM  Easting,  Northing  and  height 
of  three  widely  spaced,  non-collinear  photoidentifiable  control  points;  2)  the  x 
and  y  image  coordinates  of  the  control;  3)  focal  length  of  the  camera;  and 
4)  approximate  UTM  coordinates  of  the  exposure  station.  These  values  are 
input  into  the  calculator's  storage  registers.  The  user  then  starts  the  program 
to  compute  corrections  to  exposure  station  coordinates.  With  reasonable 
approximations  (less  than  1000  meters  away  from  the  true  values),  the  solution 
converges  after  three  or  four  iterations.  Corrections  to  the  revised  coordinates 
of  the  exposure  station  will  be  less  than  one  meter. 


The  operator  then  executes  the  command  to  compute  tilt  -  the  angle  nOp, 
imaged  nadir,  principal  point,  and  imaged  principal  point;  swing  -  the  angle  in 
the  photo  plane  measured  clockwise  about  the  optical  axis  from  the  +  y-axis  to 
the  imaged  nadir;  and  azimuth  of  the  principal  plane  -  the  UTM  azimuth  of  the 
direction  of  tilt.  These  values  are  automatically  stored  for  use  in  the  following 
algorithms. 

A  subroutine  called  Direction  Cosines  (DC)  (Appendix  B,  Line  413)  com¬ 
putes  the  vector  algebra  common  to  the  height  and  length  determinations.  The 
user  must  input  into  the  storage  registers  the  x,  y  of  the  feature  of  interest  and 
focal  length  of  the  camera.  That  is,  the  base  and  top  of  a  vertical  feature  or 
the  end  points  of  a  horizontal  feature. 

The  final  step  is  to  command  the  calculator  to  execute  the  height 
(Appendix  B,  Line  471)  or  the  length  (Appendix  B,  Line  491)  subroutine. 

The  software  was  designed  for  practical  photo  interpretation  applications 
and  will  give  results  comparable  to  those  obtained  with  conventional  vertical 
photography.  Obviously,  the  geometry  is  strongest  with  a  152.43  cm  focal 
length  camera  and  control  widely  spaced  to  form  a  strong  three-sided  pyramid. 

Distance  measurements  are  considered  to  be  small,  such  as  the  length  of  a 
bridge,  the  edge  of  a  building,  or  a  plot  boundary.  As  with  conventional  vertical 
photography,  no  compensation  is  made  for  displacement  due  to  elevation 
differences. 

5.  CONCLUSION: 

The  use  of  tilted  aerial  photography  until  now  has  been  largely  limited  to 
qualitative  applications.  A  photogrammetrist  can  use  slightly  tilted  photo¬ 
graphy  on  a  stereo  plotter  or  the  photo  interpreter  can  send  it  to  the  photo  lab 
for  rectification.  However,  these  procedures  require  elaborate  optical  equip¬ 
ment,  are  time  consuming,  and  require  highly  skilled  technicians. 


The  tools  and  procedures  for  the  digital  rectification  described  in  this 
paper  open  new  opportunities  for  the  use  of  tilted  photography.  The  compo¬ 
nents  are  inexpensive,  reliable,  and  virtually  maintenance  free.  The  system  is 
easily  learned  and  can  be  used  ip  either  the  office  or  field.  Furthermore,  a 
photo  interpreter  with  access  to  a  powerful  programmable  calculator  will 
develop  his  own  algorithms  to  make  his  time  more  productive. 
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APPENDIX  A 


Math  Model 


I.  Direction  Cosines: 

Image  coordinates  of  nadir  xn,  yn  and  zn 

x  =  pn  cos  (90  -  s) 
n 

(l) 

y  =  pn  sin  (90  -  s) 
n  r 

(2) 

HH 

1 

II 

c 

N 

(3) 

where  pn  =  f  tan  t, 
s  is  swing, 
f  is  focal  length, 
t  is  tilt,  and 

pn  is  the  vector  from  the  imaged  principal  point  to  the  imaged  nadir. 


Computation  of  direction  cosines  1,  m,  and  n. 


1,  =  Xj/dj, 


mi  •  n/di' 


"i  =  Vdi 


(4) 


where  d.  .  x.^  +  y.^  +■  (5) 

Angles  between  vectors  are  computed  using  the  cosine  law: 

<  VjOV2  =  cos~*  (lj  I2  +  mj  mj  +  n^  n2)  (6) 

where  <vjov2  *s  angle  at  the  principal  point  o  between  two  vectors. 


The  value  on  is  used  in  computing  both  heights  and  distances, 
on  =  f/cos  t  (7) 

Refer  to  Figures  1,  2,  and  3  for  graphic  representation  of  variables. 
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II.  Height  Determination: 

b*  and  t'  are  points  where  vectors  from  the  principal  point  to  the  base 
and  top  of  the  object  being  pleasured  would  intersect  a  horizontal  plane 
passing  through  the  imaged  nadir. 

nE1  =  on  tan  <  nob  (8) 

nt'  =  on  tan  <  not  (9) 


d.'  is  the  rectified  image  displacement. 

d.'  =  Tit '  -  nb'  (10) 


r'  is  the  rectified  radial  distance  from  the  imaged  nadir  to  t  ;ctified 
image  of  the  base  of  the  object. 

r'  =  nib'  (11) 


H  is  the  height  of  the  exposure  station  above  the  object. 

H  x  ZQ  -  HF  (12) 

where  Z  is  the  height  of  the  exposure  station  above  the  datum 
(computed  in  Church's  resection),  and 

Hp  is  the  elevation  of  the  object  above  the  datum. 


Substituting  the  values  into  the  conventional  height  determination 
equation 


height  s  d.'  H/r' 


(13) 


III.  Length  determination: 

The  equation  for  the  length  of  an  object  on  a  conventional  vertical 
aerial  photograph  is 

L  =  1(Zq-H)  /f  (14) 

where  L  is  the  length  of  the  object  1  is  the  imaged  length  of  the  object 
and  the  other  parameters  are  as  defined  before. 

Using  the  horizontal  plane  passing  through  the  imaged  nadir  we  can 
substitute  the  rectified  image  values: 

L  =  V  (Zq-H)  /  oH  (15) 

where  1'  is  the  rectified  image  length. 

1'  is  the  vector  subtraction  of  the  vectors  oe.*  and  o^',  where  e^'  and 
62'  are  the  rectified  image  points  of  the  ends  of  the  object  i>feing  measured. 


1' 

=  062'  -  oe^' 

(16) 

—  2  —  2 
=  oe  j  + oe2  - 

2  oe^  oe ^  cos  <e2  oej 

(17) 

oe^' 

=  on  /  sin  <  noe^ 

and 

(18) 

oe 2' 

=  on  /  sin  <  noe2 

(19) 
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FIGURE  6.  Bausch  &  bomb  Tube  Magnifier 


FIGURE  7.  Measuring  Photo  Coordinates 
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